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We review some elementary aspects of the critical properties of the series 
of metal-insulator transitions that constitute the integer quantum Hall effect. 
Numerical work has proven essential in charting out this phenomenon. Without 
being complete, we review network models that seem to capture the essentials of 
this critical phenomenon. [] 



1 To appear in Annual Reviews of Computational Physics, Vol. 5, edited by D. Stauffcr 
(World Scientific, Singapore, 1997). 
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1. Introduction 



The existence of localized states in semiconductors has been known since 
the seminal work of Anderson [1,2]. He used a deceptively simple-looking tight- 
binding model and showed that when the disorder of the binding energy on the 
scale set by the interaction energy exceeds a given limit (which is zero in one 
dimension), the wave functions are localized in space. I.e. the eigenstates fall off 
exponentially while still forming a continuous energy spectrum. The physics be- 
hind this mechanism is destructive interference. Another mechanism that leads 
to localization is that of percolation, for the first time described by Flory [3], 
and later used by Skal and Shklovskii [4] to describe hopping resistivity of lightly 
doped semiconductors. Percolation theory has, however, grown far beyond semi- 
conductor physics, and even beyond physics. For a modern review, see, e.g., 
Sahimi [5]. The review by Thouless [6] discusses percolation and Anderson local- 
ization in a unified manner. 

In this paper, we shall discuss localization in the integer quantum Hall effect. 
The classical Hall effect is standard textbook material, see e.g., Ashcroft and 
Mermin [7]. A conductor in the shape of a thin rectangular plate, i.e. a Hall 
bar, is placed in the strong perpendicular magnetic field B as shown in Fig. 1. 
A current I x flows through the bar, giving rise to a longitudinal voltage drop V x 
and a transversal voltage drop V y , which is also called the Hall voltage. It is the 
appearance of this transversal voltage that constitutes the classical Hall effect. 
The ratios 

Pxx = -j- , (1) 

and 

Py. = V f (2) 

define the longitudinal resistance and Hall resistance respectively.^] Simple argu- 
ments based on the Drude model of conductance predict the longitudinal resis- 
tance to be independent of the magnetic field B, and the Hall resistance to be 
linear in B = \B\. This is also found experimentally. 

Using a Si MOSFET inversion layer — which forces the charge carriers to 
move in a two-dimensional plane — placed in a 19 Tesla magnetic field and at 
a temperature of 1.5 K, von Klitzing et al. [8] observed a completely different 
behavior of the two resistances: The longitudinal resistance showed a series of 
pronounced peaks while being zero between them. The Hall resistance did not 
grow linearly with B, but developed a series increasingly pronounced steps. The 
reader should take a look at Paalanen et al. [9] to see these plots, which here 
are based on a GaAs heterostructure (as they did not appear explicitly in the 

■"■Note that in two dimensions, resistance and resistivity are equivalent — and that this 
system is essentially two-dimensional. 
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Figure 1: Schematic diagram of the Hall system 

original von Klitzing et al. paper). This is the integer quantum Hall effect. We 
have listed in Table 1 the typical values of some relevant physical quantities in 
the GaAs heterostructure. 

The behavior of the longitudinal resistance clearly indicates that the charge 
carriers of the quantum Hall system repeatedly localizes and delocalizes with 
increasing magnetic field. What kind of localization goes on, and why the re- 
peated localization-delocalization-localization transitions? These questions lie at 
the heart of the integer quantum Hall effect. 

Several very good texts on the integer (and fractional^) quantum Hall effect 
have already appeared, see e.g. Ref. [10-13]. Our aim with this paper is to convey 
an elementary understanding of the phenomenon, rather than a full-scale review 
of the entire literature. We therefore urge the readers to consult these more 
complete works. 

2 Which appears in clean samples at even higher magnetic fields and lower temperatures than 
are necessary for the integer effect, in the form of steps appearing with a different regularity 
than those of the integer quantum Hall effect. 
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We should perhaps apologize that the parts of this review that deal with the 
numerical studies of the integer quantum Hall effect is short compared to those 
dealing with describing the physics of the system. However, we believe this to 
be justified. The numerical models are simple, but the results that they lead to 
are profound. In order for the reader to appreciate the numerical models fully, 
we have chosen to emphasize the physics of the effect. We write this review for 
an audience in mind with a background in statistical 



m = 0.067m e 


Effective mass 


n = 4- 10 15 m" 2 = (16nm)- 2 


Electron sheet density 


Ep = 14meV 


Fermi energy 


Xf — 40nm 


Fermi wavelength 


A = 10 2 - 10 4 nm 


Elastic mean free path 


A^ > A 


Phase coherence length 


I = ^h/eB = 26( J B/T)- 1 / 2 nm 


Magnetic length 



Table 1: Some relevant physical quantities in the GaAs heterostructure 



physics. We have therefore attempted to be quite thorough and explicit in 
setting up the physical framework. No prior knowledge to the problem of two- 
dimensional quantum transport is assumed, and only elementary quantum me- 
chanics is used. The essence of what will follow is this: The integer quantum 
Hall effect is closely related to the Landau levels of the electrons, formed when 
the system is placed in a strong perpendicular magnetic field. The energy of the 
electrons is quantized in way similar to the harmonic oscillator, see Eqs. (13) and 
(15). Whenever the magnetic field is reduced and one of these Landau levels in 
the bulk sinks through the Fermi level of the electrons, an additional contribution 
to the electrical conductance is made possible by filling this new Landau level 
with electrons from the attached reservoirs. This leads to the quantized current 
of Eq. (68) and to the subsequent "hydrodynamical" scenario thereafter in which 
disorder plays a crucial role. 

The organization of the review is as follows: We discuss in the next section 
the behavior of a two-dimensional free electron gas in a strong perpendicular 
magnetic field. The integer quantum Hall effect enters in a very subtle way 
in this simple model. When we introduce interactions between the underlying 
lattice and the electrons, the integer quantum Hall effect becomes much more 
robust — and easy to understand. We discuss this in Sec. 3. As will be evident 
from this section, percolation plays an important role. In Sec. 4, we discuss 
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the critical behavior of the integer quantum Hall system, which is caused by 
localization-delocalization transitions when the the externally applied magnetic 
field is changed. We furthermore discuss a famous semi-classical argument that 
seemingly identifies the nature and thus the universality class of these transitions. 
Unfortunately, there are problems with this argument, as we point out. In Sec. 
5, we discuss a class of numerical network models that seems to capture the 
critical behavior of the quantum Hall system very well. Section 6. contains our 
conclusions. 
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2. The Free Electron Gas in a Closed System 

The Hamiltonian of an electron moving in the (x, y) plane and in a perpen- 
dicular, constant magnetic field B = Be z is 

H = ±(p+eA)\ (3) 

— * — * 

where B = V x i, e is the elementary charge and m the effective mass. We 
shall temporarily not specify any particular gauge. The standard commutation 
relations apply: 

[x,Px] = [y,Py\ = ih , (4) 

and 

[x,y] = [Px,p y ] = [x,p y ] = [y,p x ] = . (5) 

Following Kubo et al. [14], we perform a canonical transformation from the orig- 
inal coordinates (x,y) and (p x ,P y ) to the relative coordinates (£,77) and guiding 
center coordinates (X,Y): 

i = ^{p y + eA y ) , (6) 

V = (Px + eA x ) , (7) 
X = x-i, (8) 

and 

Y = y-r 1 . (9) 
The commutation relations (4) and (5) lead to 

[X,Y] = -[^, V ]= i (^ =il\ (10) 

and 

% X] = [77, Y] = K, Y] = [77, X] = . (11) 
In Eq. (10), we have used the magnetic length 

12 = h ■ (12) 

which corresponds to the cyclotron radius for the lowest quantum state. The 
accompanying classical cyclotron frequency is 

u= 6 -^. (13) 
m 
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The relative coordinates £ and r] describe the position of the electron relative to 
the center of the circle on which the electron is moving. The position of this 
center, the guiding center, is given by X and Y. The Hamiltonian is in these 
coordinates 

Ml 

(14) 



m 



which we recognize as the harmonic oscillator Hamiltonian by appealing to the 
commutation relation (10). The energy spectrum is therefore 



E„ = Uuj I n 



1 



(15) 



These are the Landau levels. 

In order to construct the wave functions, we must specify a gauge. We choose 
the Landau gauge where A x = (—By, 0). The Schrodinger equation then becomes 



2m 



dy 2 



' d ieB N 
dx h ^ ' . 



i) x n = E n i) x n , 



(16) 



where the superscript x refers to the choice of gauge that we have made. We now 
set 



so that Eq. (16) becomes 

h 2 d 2 



m , 

- H uo z 

2m dy 2 2 



Xn{y) , 



Xn(y) = Kx(y) 



:i7) 



(18) 



Thus, the energy eigenfunctions ip x Y consist of plane waves in the positive x 
direction and harmonic oscillator wave functions Xn(y) — Xn(y — kl 2 ) in the y 
direction centered around Y = kl 2 . The energy levels are independent of k. The 
energy eigenfunctions are therefore 



e iYx/l 2 



(19) 



These eigenfunctions form a complete set, and any wave function may be written 
r(x, y) = E / dY^L X n(y - Y)C n (Y) . (20) 



A change of gauge A — > A— Va leads to a change in phase of the wave function, 
■ip — > exp(iea/h)ip. For example, we may go from Landau gauge A x = (—By, 0) 
to the symmetric gauge A s = (—By/2, Bx/2). Integrating the difference between 
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these, Va = (—By/2,—Bx/2), gives a = —Bxy/2. An arbitrary wave function 
in this gauge may therefore be written 

. -i(xy-2Yx)/2l 2 

r(x, */) = £/ dY -= X n(y - Y)C n (Y) . (21) 

n J V L 

Let us now count the number of states per Landau level and per area in the 
system. Assuming periodic boundary conditions in the x direction, the wave 
function (17) must be periodic in the length L of the system, ^ Y {x + L,y) = 
i/j* y (x, y). This quantizes the possible values of Y, 

2nl 2 . . 

Y = —3> ( 22 ) 

where j = 0, ±1,±2, The system has width W in the y direction. Thus, 
< Y < W — remember that Y determines the center of the harmonic oscillator 
wave functions in the y direction.^ Combination of this restriction with Eq. (22), 
leads to j being confined to the integers on the interval < j < LW/2tt1 2 = j max . 
Thus, the number of states per Landau level per area is 

.7 max 1 / n o \ 

nB= LW = 2^- (23) 

We now turn on an electric field £ in the y direction. The Hamiltonian (14) 
then becomes 

H = juj 2 (e + V 2 ) + e£ y (Y + V ) , (24) 
where we have used Eq. (9), y — Y + rj. We may rewrite this Hamiltonian as 

2c2 



From this it is possible to write down the eigenf unctions and corresponding energy 
levels directly by inspection: 

rn,y(^y) = -^Xn(y-Y + ^L) , (26) 

and 

^^.( n + i) + ^(y-^) + f(|) 2 , (27) 

where the second term represents potential energy in the electric field and the 
third term the kinetic energy associated with guiding-center drift. 



3 We ask for patience from those who feel that our treatment of the boundaries of the system 
in the y direction to be too simple minded. We will return to these boundaries in great detail 
in Sec. 3.2. 
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Let us now calculate the current densities in the system.^ Essential ingredients 
of such a calculation are the expectation values of the velocity in the x and y 
directions, 

v x (n,Y) = dx dy ip x n Y *(x, y)xip x n Y (x, y) , (28) 



and 

v v (n, Y) = J dx dy < >y *(x, y)yip* >Y (x, y) . (29) 

We split x and y into x = £ + X and y = rj + Y . The Heisenberg equations of 
motion give 

d£ i [TT M ( e£ y 



and 



* n lH ^ = -^" + ^)' (30) 

*- = ip,,]=4«{, (31) 
dX i el 2 £ 

f = = T^ = §- (32) 

f = i[H,y]=0. (33) 

We have here used the commutation relations Eqs. (10) and (11), and the Hamil- 
tonian, Eq. (24). Combining these results with equations (28) and (29) give 

v x (n,Y) = $L, (34) 

and 

v v (n,Y) = 0. (35) 

In Eq. (34) we have used that 



dr] \Xn(v + eS y /muj 2 )\ 2 (r] + e£ y /mu 2 ) = / drj \Xn(v)^V = • (36) 

3 

In equation (35), we have used that 



dv Xn(v + e£y/mu )*^Xn(v + e£ y /mu ) = 

3 

dr) Xn{v + e£y/muj 2 )*(-il 2 )--^Xn(v + e£ y /muj 2 ) = 0. (37) 

dr] 



4 By using the expression 

1 dE n , Y 

v ~ n dk ' 

where E nY is given by Eq. (27), we obtain the results Eqs. (34) and (35) immediately. However, 
the way we present the derivation of Eqs. (34) and (35) is more detailed, and is in fact a 
derivation of the above equation. 
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With n e electrons per area, Eqs. (34) and (35) give rise to a current density 
in the x direction equal to 

e 2 n e 2 e^"u 
jx = - ev A n i Y ) = -2irl 2 n e — £ y = - — £ y = — — £ y . (38) 

The sum runs over all electrons and we ignore spin. In the y direction the current 
density is 

jy = -T, ev v^ Y )=°- ( 39 ) 

n,Y 

We have in Eq. (38) defined the filling factor 

v = ^. (40) 
n B 

lis is defined in Eq. (23). 

From Eqs. (38) and (39), we read off the elements of the conductivity 

0~XX = ^yy = , (41) 

and 

e 2 

0~ yx = -&xy = -J^ v ■ (42) 

From this the resistance tensor follows by inversion, or directly by reference to 
Eqs. (1) or (2) as 

Pxx = Pyy = ® ■ (43) 

and 

Pxy = ~Pyx = ~~2^j ' (44) 

The resistances (43) and (44) capture one aspect of what is measured in the 
integer quantum Hall effect: zero longitudinal resistance. However, the steps, 
which have their centers at integer values of the filling factor v and the onset 
of longitudinal conductance when the system moves from one step to the next, 
are not found in this system. The reason that the steps are not visible, is that 
the system is closed. In Sec. 3.2, we will consider an open system, where steps 
do appear. The reason why steps only appear in open systems will be clear at 
that point. Before we turn to this discussion, we add interactions between the 
electrons and the underlying system, as they play a crucial role in what happens 
between the steps — which is the subject of this review. 

3. Interactions between Electrons and Lattice: The Disordered Potential 

We now add an arbitrary potential V(x,y) to the Hamiltonian (3). The 
Schrodinger equation we have to solve is then 

Hij)(x, y) = H ip(x, y) + V(x, y)^(x, y) = Eip(x, y) . (45) 
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We expand the wave function as in Eq. (20) In this basis, we have that 



H r(x,y) = (n> + ±) J dY' ^ Y ' xjl \ nl {y - Y')C n ,(Y>)/VL . (46) 

By multiplying Eq. (45) by e~ tYx ^ 2 Xn(y — Y) / ^/Z and integrating over (x,y), the 
Schrodinger equation is transformed into 

J dx dy V(x,y)Y: J dY' e^'-W^y - Y) Xn >(y ~ Y')C n ,(Y')/27rl 2 

n' 

— {E — E n )C n (Y) . (47) 

This expression is similar, but not identical, to the one given by Tsukada [15], 
which is the standard reference, but contains a small error. 

Let us now make some changes of variables in Eq. (47). By defining Ay = 
Y' — Y and using that x = X + £ and y = Y + r], we can write Eq. (47) as 

- , iA Y (X+0/l 2 

J d£d v V(X + Z,Y + V )J dA Y — 

J2xUv)Xniv-^y)CniY + A Y ) = (E — E n )C n (Y) . (48) 

n' 

So far, we have made no approximations. However, from now on we shall 
assume that V(x, y) is slowly varying on the scale of the magnetic length /. Since 
the functions Xn(x) are localized over distances of the order of / — they are 
harmonic oscillator eigenfunctions — we may substitute V(x,y) by V(X,Y) in 
Eq. (48) as the lowest order term in an expansion in I. With this substitution, 
Eq. (48) simplifies dramatically by first integrating over £ and then using the 
orthonormality of the basis functions Xn, 

V(X, Y)C n {Y) = (E- E n )C n (Y) . (49) 

The corresponding energy eigenfunctions are 

,. iYx/l 2 

r n (x, y) = J dY—j=-Xn{y - Y)C n {Y) . (50) 

What is the physical contents of Eq. (49)? In the classical limit, where X and 
Y are treated as c-numbers, it predicts a motion that follows the equipotential 
curves of V(x, y) = E — E n . The classical electrons move in a strong magnetic 
field such that they neither gain nor loose energy. In order to expand on this, we 
write down the Heisenberg equations of motion for X and Y: 

dX i ITT VI 1 dV 

w-fF'^+rtdj;' (51) 
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and 

ir = h [H > Y] = -^B^> (52) 

where H = Hq + V, where Hq is given in Eq. (3), and we have used the commu- 
tation relations (10) and (11).Q There are no approximations in Eqs. (51) and 
(52). Introducing the same level of approximation as was done to arrive at Eq. 
(49), i.e., assuming that V is slowly varying on the scale of the magnetic length 
/, Eqs. (51) and (52) become^] 

dX 1 dV 

^ = + 7bW> (53) 

and 

dY^ J_ dV_ 

dt ~ eB 8X ' { ' 

These equations are also valid in the classical limit. Let us now multiply the 
velocity vector of the guiding centers (X, Y) with the gradient of the potential, 

W = (dV/dX, OV/dY). We then get U ■ W = by use of Eqs. (53) and (54). 
Here U = (X, Y). Thus, the motion of the guiding center is always orthogonal to 
the local gradient of the interaction potential, and as a consequence, the guiding 
centers follow the equipotential curves of V. 

3.1 Motion Near a Saddle Point: Tunneling 

There will be quantum mechanical tunneling if two distinct equipotential lines 
come close at some point in space, as shown in Fig. 2. This may be somewhat 
surprising in that the tunneling has to appear in the direction orthogonal to the 
paths. We may understand this qualitatively through the following semi-classical 
picture: The paths drawn as in Fig. 2 are those of the guiding centers, which 
the electron itself is rotating about. Thus, the momentum vector of the electron 
points in "all" directions, even orthogonal to the path that the guiding center is 
following — making tunneling possible in the seemingly most unlikely direction. 

Let us analyse this transverse tunneling more closely. Following Fertig [19] 
and Mil'nikov and Sokolov [20], we concentrate on the region near a saddle point 
of V. We choose a coordinate system so that the saddle point is placed at the 



5 Write V as a Taylor series in x and y, then calculate the commutators term by term and 
reassemble the result. 

6 If the potential V is not caused by the surrounding medium in which the electron move, 
but rather their mutual interactions, Eqs. (53) and (54) become in the classical limit (when X 
and Y are c-numbers) the Kirchhoff equations of motion for vortices. These have many strange 
properties, for example that the three-body problem is integrable, see Grobli [16], Synge [17] 
and Aref [18]. 



12 



Figure 2: The potential V has a saddle point in region A. Even if the two paths 
marked in bold, which follow the equipotential curves of V, are antiparallel where 
they are at the closest, there is tunneling between them. 

origin. Furthermore, we orient the coordinate system so that close to the origin, 
we have that 

V(X,Y) = V -aX 2 + f3Y 2 , (55) 

where a > and (3 > 0. Thus, V has a local maximum at the origin when moving 
along the X axis, and a local minimum at the origin when moving along the Y 
axis. Two situations may then occur: (1) E — E n < V(0,0) = Vq, where E is 
found by solving the Schrodinger equation (49). The corresponding semi-classical 
paths are shown in Fig. 3a. Likewise, the situation when E — E n > Vq is shown 
in Fig. 3b. Suppose that the situation is that of Fig. 3a. Then, it is convenient to 
use the representation Y = il 2 d/dX, thus satisfying the commutator [X, Y] = il 2 , 
Eq. (10). Equation (49) then becomes^] 

-/3/ 4 ^ - al 2 j C n {X) = (E-E n - V ) C n (X) . (56) 

If, on the other hand, the situation is the one described in Fig. 3b, we switch 
to the representation X = —il 2 d/dY, which also satisfies the commutation rule 
[X, Y] = il 2 . Equation (49) then takes the form, after an overall change of sign, 

-«/ 4 ^ - PY 2 ^ C n (Y) = (V + E n — E) C n (Y) . (57) 

Thus, we are faced with solving the one-dimensional scattering problems posed 
by the Schrodinger Eqs. (56) and (57) in order to study tunneling. This is difficult 
since \V\ — > oo when \X, Y\ — > oo. However, we can get around this difficulty 
by substituting V - aX 2 -> V / cosh 2 (^a/V X) in Eq. (56) or -V Q + (3Y 2 -> 
Vq/ cosh 2 (^J j3/V Y) in Eq. (57), see Hansen [21]. This substitution leads to the 

7 Since the commutator is of O(P), one could worry about the consistency of Eq. (56). A 
tedious check leads to a reassuring result: The only extra term needed for consistency represents 
a constant shift in energy near the saddle point. We ignore this irrelevant constant in what 
follows. 
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a 



b 



Figure 3: The potential V has a saddle point at the origin. In (a) we show two 
semi-classical paths corresponding to E < E n + V , and in (b) we show the paths 
corresponding to E > E n + V . 



correct quadratic behavior near the origin, and plane-wave solutions far away 
from the origin. The exact solution of the resulting one- dimensional scattering 
problem may be found, e.g., in ter Haar [22] or Landau and Lifshitz [23]. Both 
Eqs. (56) and (57) lead to a transmission coefficient 



T = , (58) 

1 + P K ' 



where 



P = e -27r/l 2 y/T/^\E-E n -V \ _ ( 5 g) 

We make a couple of comments: (1) When AE n = (E — E n — Vq) — > 0, the 
transmission coefficient T approaches the value 1/2, which is the maximum value 
it may take. Thus, when AE n = 0, an electron approaching the saddle point 
will leave it following either of the two possible paths with equal probability. (2) 
Using the WKB approximation, as was the approach in Mil'nikov and Sokolov 
[20], leads to T WKB = P = T/(l - T) = T + T 2 + • • •. 

The potential V(x, y) is caused by disorder in or near the inversion layer of the 
MOSFET, which in turn may be induced by, e.g., crystal defects or impurities. 
Thus, the potential may be seen as a random landscape of valleys and ridges. 
Electrons move along equipotential curves in this landscape. Most of these curves 
form closed loops, which may join across saddle points. How can such a scenario 
give rise to the strange behavior of the conductivity that one observes in the 
integer quantum Hall effect? Before attempting to answer this question, we 
discuss the concept of conductivity in the present context. 

3.2 Conductivity and Edge States 

Conductivity is a local quantity. Conductance is not. In the following discus- 
sion, it will appear that conductance is the fundamental quantity in the quantum 
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Hall system, rather than the conductivity. However, this may be misleading, and 
just a result of the framework we use to present the results. To make an analogy, 
presenting the Maxwell equations in the form of surface and line integrals makes 
them look very different, but the physics is of course the same. We follow in this 
section mainly the work of Biittiker [24]. 

So far, we have not cared about the edges in our Hall system of size L x W. 
The only place where they have entered the discussion, is where we counted the 
number of states per area. Let us look at the solutions of the Schrodinger Eq. 
(45) for the Hall system confined to a bar of finite size, see Halperin [25]. First, 
we assume that there is no disorder potential V. We have already determined 
the energy eigenfunctions for the infinite system, Eq. (19). These eigenfunctions, 
if)* Y (x, y), consist of "stripes" oriented in the x direction of width ~ I. The 
position of the middle of a given "stripe" is determined by the quantum number 
Y . The wave functions must have nodes at y = and y = W, i.e., at the edges of 
the Hall bar. Thus, for Y = or Y = W, only the odd harmonic oscillator wave 
functions, X2m+i(y) fulfil the boundary condition. Thus, 2m + l = l=^m = 
corresponds to the lowest energy eigenstate here, and the corresponding energy 
level is Eq y=o — 3hu/2. Moving away from the edge and into the sample, i.e., 
increasing Y, the energy levels move continuously to their bulk values, which 
for the ground state is Eo,y — hu/2. This is general: The energy eigenvalue 
corresponding to the quantum number n takes on the value E n y = hu(n + 1/2) 
inside the sample, but rises sharply to the values E n> Y=o = E n y =w = fuo{2n + 
3/2) at the edges. This is illustrated in Fig. 4. It should be noted in this figure 
that the energy levels E n y is continued for Y values that are less than or 
larger than W. This makes sense. Classically, this means only that the center 
of the semicircles that the electrons follow lie outside the sample. The charges 
themselves, whose coordinates are (x,y), never move outside the sample limited 
by y = and y = W. 

Are there any currents flowing in this system? Judging from Eqs. (38) and 
(39), the answer is no when there is no external electric field S y present. However, 
this is not true, as we did not analyse what is happening at the edges of the 
sample. From Eqs. (28) and (36), we have thatQ 

r w 

v x (n,Y) = -co dy\x n . Y (y-Y)\ 2 (y-Y) , (60) 
J o 

where X n>Y (y) X2n+i(y) when Y -> or W, and Xn,v{y) -> Xn{y) when Y > 
and Y <C W. The wave function Xn falls off over a range I on each side of Y, 
except at the borders, where it is zero outside the Hall bar. Thus, when F>0 
and Y W, Eq. (60) integrates to zero — as was already used in Eq. (36). 

8 We are implicitly assuming periodic boundary conditions in the x direction by our choice 
of wave functions, Eq. (19). Thus, we only discuss the boundaries y = and y = W . This is 
known in the literature as Corbino ring geometry. 
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However, when Y — 0, Eq. (60) leads to 

v x (n,0) = -uj dy \x2n+i(y)\ 2 y < , (61) 
Jo 

and for Y = W, 

rW 

v x (n, W) = -u dy \ X 2n + i(y - W)\ 2 (y - W) = -v x (n, 0) . (62) 



There are therefore currents flowing in opposite directions along the edges of the 
system. 

Let us now attempt to rederive the conductances, Eqs. (41) and (42), but 
without applying an external field £ y explicitly. The reason for this is that it is 
awkward to work with small but finite external fields when the disorder potential 
V(x,y) is present. And the reason why one should believe that to be possible, is 
that conductance is an equilibrium property of the system. This will be clear in 
a moment. We shall find that we fail slightly, but in an interesting way. 

Let us calculate 

-jy- = l dyx n ,y(y- Y )^xnAy-Y) 
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f w dH 

J o d y x*n,v{y - Y )^yXnAv- Y ) 

[dy r n , Y (y - y)^ I? + & ~ Y » 2 } ^ y - y ) 



rW 

-moj 2 



/ dyx*nAy- Y )[y- Y ]xnAy-y)- (63) 

Jo 

Finally, by combining this expression with Eq. (61), we derive the dispersion 
relation 

^, y) j;^. (64) 

Once again, we see that currents only occur at the edges, which since only here 
E n: Y is an explicit function of Y, see Fig. 4. 

We now open the system by cutting the periodic boundaries in the x direction 
and connect the Hall bar to two reservoirs. One is kept at a chemical potential 
/ia, and the other at a chemical potential hb- This is shown in figure Fig. 5. 
What is the total current flowing from one reservoir to the other one through 
the Hall bar? We start by only considering the current associated with the nth 
Landau level, I XjU . This current is 

/■Ms 

I n ,x= dp ev x {n,Y(n))g n 1ji) , (65) 

J )1A 

where we use the one-dimensional density of states propagating to the right, 

1 d Jn (E) 1 dY n (E) 
9n{E) = L^E-=2^^E-- (66) 

The reason for using the one- dimensional expression, is that the edge currents 
are one-dimensional. We have used Eq. (22) in going from j to Y on the right 
hand side of Eq. (66). The limits of the integral in Eq. (65) reflect that it is 
the excess current we are measuring, as the currents associated with matching 
chemical potentials cancel (but they flow along opposite edges). Combining these 
two last equations with Eq. (64), we find that the total current in the nth Landau 
level is 

I n ,x = -^{fiB- Ha) ■ (67) 

These currents are, as we have seen, edge currents. There is one such current for 
every bulk Landau level which is below the Fermi level. The reason for this may 
be seen in Fig. 4: The Landau levels bend upwards at the edges, and those below 
the Fermi level in the bulk cut through the Fermi level at the edges. Thus, the 
total current is, when both (j,a and /is differ infinitesimally from the Fermi level, 

I x = N^(fji B - fi A ) , (68) 

where N is the number of bulk Landau levels below the Fermi level. It is an 
integer. 
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The factor Ne/h in Eq. (68) is an inverse contact resistance (divided by e), 
meaning that it is due to a redistribution of the population of states when the 
charges enter the system from that they had in the reservoirs. If one measures 
the voltage difference between contacts 1 and 2, or 3 and 4, which are shown in 
Fig. 5, one finds zero: There is no change in chemical potential along the edges. 
However, measuring voltage difference across the system, i.e., between contacts 
1 and 3, 1 and 4, or 2 and 3, or 2 and 4, one finds a potential difference fiB — Ha- 
Thus, there is no resistance in the flow direction of the current, and there is a 
resistance in the direction perpendicular to the flow direction. In other words, 
we have 

Pxx = Pyy = ® , (69) 

and0 

h 

Pyx = -Pxy = ■ (70) 

These resistances should be compared to those listed in Eqs. (43) and (44). The 
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Figure 5: A clean Hall bar (i.e., V = 0) in contact with two reservoirs A and 
B kept at chemical potentials ^a and Furthermore, four contacts have been 
added so that the potential between any pair among them may be measured. 
There is no current passing through contacts 1, 2, 3, and 4. 

difference between the two sets of resistances occurs in Eqs. (70) and (44): In Eq. 
(70), the Hall resistance is quantized, while in Eq. (44), it is not. 

Why this difference? In Sec. 2, we worked with a closed system. The number 
of electrons was fixed. The present system is open; it is in contact with two 

9 The extra 1/e factor in Eq. (70) comes from measuring chemical potential — which has 
the units of energy — in units of electric potential. 
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reservoirs. Whenever a new Landau level drops below the Fermi level, e.g. through 
decrease of the magnetic field B, it fills up due to the contact with the reservoirs. 
This is not possible in the closed system. 

We now "turn" on the disordered potential V. The arguments leading to the 
bending upwards of the Landau levels at the edges of the clean Hall bar rested 
on there being sharp boundaries on the scale of the magnetic length I, cutting 
the wave function off. It is more realistic to assume that the boundaries are 
sharp on the scale of the fluctuations of the disorder potential V, which in strong 
magnetic fields is much larger than the magnetic length /. Thus, at the edges, 
the potential V turns smoothly upwards on the scale of /, while the Landau levels 
remain unchanged throughout the system.^ 

We have seen in Sec. 3 that the electrons follow equipotential curves in the 
energy landscape. In the bulk of the Hall bar, most such curves form closed loops 
circling "mountains" or "valleys" in the energy landscape. However, at the edges, 
there are equipotential paths reaching from one end of the sample to the other 
due to the upwards bend of the energy eigenvalues. Connecting the Hall bar with 
reservoirs A and B kept at chemical potentials \xa and \is will induce currents 
along these paths. How many such paths are there? The argument is exactly the 
same as in the case with V = (i.e., the clean case): There is one path for every 
Landau level that is below the Fermi level in bulk. 

It is useful to picture the following situation: Take a plate and with a hammer 
knock it into a shape that resembles our disorder potential. At the edges, bend 
the plate upwards. This correponds to the bending at the edges of the Landau 
levels. Then make a series of copies of this plate and stack them at regular height 
intervals. This is a three-dimensional model of the energy levels in the Hall bar. 
Now, place the whole construction in a tub and start filling it with water. The 
water level represents the Fermi level. As the water rises, "shore lines" will form 
where the dry plates dives into the water. These shore lines correspond to the 
equipotential curves at the Fermi level. We see that for most water levels, shore 
lines run parallel to upwardly bent parts of the plates but not across the plates, 
orthogonal to the bends. However, for certain levels, there are connected shore 
lines running across the sample in connecting the two bent edges. This happens 
every time a new plate (i.e., Landau level) sinks into the water (i.e., the Fermi 
sea). 

Returning from the tub to the Hall bar, how does this appearance of paths 
running across the sample parallel to the y direction (i.e., from edge to edge) 
affect the conductances of the system? Before we answer this question, we must 
determine the conductances of the system as long as such paths are not present. 
In order to argue this, we note that the arguments leading to Eq. (68) are more 
general than the derivation seems to indicate. The point is that this conductance 

10 The conclusions we draw in the following are valid even when the edges are sharp on the 
scale of I — see Halperin [25] . 
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is, as already pointed out, a contact conductance caused by the redistribution of 
the Fermi distribution in the reservoirs at the contacts between the reservoirs and 
the conductor. Biittiker [24] therefore simply assumes that each current channel 
(i.e., edge state) is a one-dimensional perfect conductor. The Fermi energy is 
in the conductor Ep = E n + h 2 k 2 /2m, where m is the effective mass of the 
electrons. The velocity of the electrons is v n = dE n /dk. The density of states 
g n (E) = dNn/dE = (l/2ii)dk/dE, since dj\f n /dk = 1/27T in one dimension. Thus, 
the combination ev n (E)g n (E) = e(dE/dk)(dk/dE)(l/27ih) = e/h. The current 
in channel n is I n = dfiev n (fi)g n (fi) = (e//i)(/ie — Ma) 3 i- e -> Eq. (68). As 
a result, we find longitudinal and transversal (i.e., Hall) resistances to be those 
listed in Eqs. (69) and (70). □ 

In Fig. 6, we show a "dirty" Hall bar connected to two reservoirs A and B, and 
four contacts 1, 2, 3, and 4. Following the Landauer-Biittiker formalism for phase 
coherent transport [26], we write down the balance equations for the 6-terminal 
system of Fig. 6. There is a current associated with contact i in this system. 
The sign is chosen so that it is positive if it enters the Hall bar. We measure the 
chemical potentials at the reservoirs A and B in units of voltage, fi A = gV a and 
/Ab — eV#. We assume that the Fermi level is set so that iV Landau levels are 
beneath it, while Landau level number N + 1 has equipotential curves cutting 
across it from edge to edge at the Fermi energy. There is a transmission coefficient 
T associated with reaching contact 1 from contact 2 and reaching contact 4 from 
contact 3, and a reflection coefficient R associated with reaching contact 4 from 
contact 2 and with reaching contact 1 from contact 3. Furthermore, R + T = 1. 
Due to the symmetry of the problem under a change of direction of the magnetic 
field, B — > —B, there is only one reflection and transmission coefficient, and not 
two as in the more general case. We write down the balance equations, 



and 



I A = 


j((N + l)V A 


~(N+l)Vx) 


•) 


(71) 


Ib = 


j((N + l)V B 


-(N + 1)V A ) 


1 


(72) 


h = 


j((N + W- 


- (N + T)V 2 - 


- RV 3 ) , 


(73) 


h = 




-(N + 1)V B ) 


) 


(74) 


h = 




-(N + 1)V A ) 


) 


(75) 


h = 




(N + T)V 3 - 


RV 2 ) ■ 


(76) 



n Biittiker [24] discusses inelastic scattering, and shows that even in this case, the resistivities 
will be given by Eqs. (69) and (70) as long as there are no scattering paths connecting the edges 
across the sample. 
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The experimental configuration that allows one to measure the longitudinal and 
Hall resistances is one where 

h = h = h = h = . (77) 

It immediately follows that 

I A = -I B = I . (78) 

We now choose one of the six potentials equal to zero, e.g. Vb = 0, to anchor the 
voltage scale. Equations (71) - (76) may then be transformed into 

/ = j(N + T)V A , (79) 

V 2 = 0, (81) 

V 3 = jV A , (82) 

and 

e 2 N + T 

v > - TnTT v * ■ < 83 > 

We may now determine the resistances between all pairs of contacts: 

R ^-^r~ ~ ^(n+i)(n+t) ' (84) 
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R ^~— — - 7 2 (N+1)(N + T) ' (85) 
V4 — V2 h 1 

^ = ^ = , (87) 



■41 

and 

^ 32 " - F" " e 2 (TVTT) • (89) 

If now T = 1 and i? = 0, which is the case when there is no path across the 
sample at the Fermi level, Eqs. (84) to (89) reduce to 

Rl2 = #34 = , (90) 

and 

i? 42 = i? 3 i = -R41 = -R32 = -r . (91) 

e z T\ + 1 

This is the integer quantum Hall effect, and we recognize Eqs. (69) and (70) in 
the last two equations. 

From Eqs. (84) to (89) we also see what happens when the system moves 
from one Hall plateau, i.e., when Eqs. (90) and (91) reign, to the next: A longi- 
tudinal resistance appears, as long as R = 1 — T > 0. This is the localization- 
delocalization transition that appears in the Introduction to this article. And, 
we see now what causes it, namely the appearance of paths across the sample in 
the y direction at the Fermi level. But, what kind of transition is this? In the 
early eighties it was suggested that it is the same mechanism as that of coalescing 
ponds during rainy weather, i.e., it is a percolation transition, see Kazarinov and 
Luryi [27] and Trugman [28]. However, as has now become clear, this picture was 
too simple minded. 

4. The Critical Points 



After having set the scene, we are now ready to discuss the phenomenon which 
is the aspect of the integer quantum Hall effect that forms the focus of this review: 
The critical behavior associated with changing the quantized filling factor from 
N to 7V± 1. 

4.1 Expectations and Observations 

As stated at the end of Sec. 3, it was early theorized that the critical behavior 
is that of "standard" percolation. The theory goes as follows: 
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There are three classes of electron wave functions in the system: (i) Wave 
functions localized to closed paths that follow the sides of "valleys" in the energy 
landscape V, (ii) those localized to closed paths circling the "mountains", and 
(hi) those localized to the group of paths that consists of paths of type (i) or (ii) 
and which are connected through saddle points. For one particular level V(x) = 
constant, there is at least one connected cluster of paths (i.e., a separatrix of type 
(hi)) that spans the system, connecting opposite edge states. The corresponding 
energy level we denote E c ^ n . As explained in Sec. 3.2, a longitudinal conductance 
appears each time the Fermi level passes through an E c ^ n by for example adjusting 
the perpendicular magnetic field B. There is a particular magnetic field B C}71 
corresponding to the energies E c>n . 

The typical (linear) size of area the connected wave functions span, £ p , diverges 
asQ 

£ p <x\E-E c>n \- v <x\B-B c , n \~ v , (92) 

where the percolation connectivity length exponent v is known to be 4/3 in two 
dimensions. The localization length, £, is the length scale over which the wave 
function does not fall off faster than algebraically, and behaves as 

foe \E-E Ctn \-» , (93) 

where v is the localization length exponent. If percolation, in the sense described 
above, were the only mechanism, we would have 

4 

v = v = - , wrong! (94) 
3 

This, however, is not what is observed in experiments, nor in computer simula- 
tions. Koch et al. [29,30] measured v = 2.3 ±0.1, while Wei al. [31,32] found the 
value v = 2.4 ± 0.2. This value was not measured directly, but relies on assumed 
values of other exponents. On the other hand, Shashkin et al. [33,34] and Dol- 
gopolov et al. [35] observe which is not too far from the percolation value 

v = 4/3. However, also in these three latter papers, v was not measured directly, 
but relied on theoretical arguments. 

Several numerical simulations have been performed to determine v. Models 
based on continuum descriptions of the problem, using long-wavelength Gaussian 
potentials or short-range impurity potentials, result in v « 2.3 - 2.4 [36-41]. 

Network models, originally introduced by Chalker and Coddington [42] give 
a similar value for u, see [43,44]. We describe these network models in detail in 
Sec. 5. 

These for the most part mutually consistent results strongly hint at the 
localization-delocalization transition not being describable by classical percola- 
tion alone. Quantum effects play a role. 

12 The last proportionality in this equation is true as long as E = E{B) is an analytic function 
near -B Cj „, defined as E c _ n — E(B c ^ n ). 
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There are two ways that quantum mechanics may change the universality 
class of the transition, while still keeping the physical picture developed in the 
preceeding sections of this review and leading to the percolation picture just 
presented: 1) Tunneling at the saddle points and 2) interference at the saddle 
points, making this an example of Anderson localization. 

4.2 The Mil'nikov and Sokolov Argument 

There is a famous argument, first presented by Mil'nikov and Sokolov [20], 
that leads to 

u = u+l = -, (95) 

based on purely semi-classical reasoning, i.e., taking only tunneling into account 
while ignoring interference. 

We follow first Hansen and Liitken [45] in their interpretation of the Mil'nikov- 
Sokolov argument. A different approach, which we also sketch, can be found in 
Zhao and Feng [46]. 

Let us assume that there is a typical distance a between saddle points. We 
tune the magnetic field, and find classical overlap at saddle points in a window 
AB around B cn which creates a connected cluster across the sample. Tunneling 
changes the question of overlap: Overlap is then no longer characterized as overlap 
within the magnetic length /, but by the tunneling probability T given by Eqs. 
(58) and (59). Let us write Eq. (59) as 



where 



P = e~ 1/x , (96) 



X he\AB\ ■ 1 ' 



Classically, we may define an "elementary" localization length as the typical 
distance between saddle points, a. This is the length over which the electrons 
may move if there are no overlaps at the saddle points, as they are constrained 
to follow closed paths with this average diameter. 

When tunneling is included, the elementary localization length is changed. 
When AB is small, we have from Eqs. (58) and (97) that T ps e" 1 ^. An electron 
starting somewhere in the system will pass K such saddle points with a proba- 
bility e~ K / x (when interference effects are ignored). Fixing this probability to, 
say, e^ 1 , we may define a tunnel length to be 

A = Ka = x a ■ (98) 

From Eq. (97) we see that A ~ |Ai?| _1 when B — > B c . In terms of the percolation 
picture, A takes over from a as the lower cutoff in length scale when tunneling 
is taken into account, an assumption relying on the one-dimensional character of 
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the percolation clusters near threshold.0 The percolation correlation length £ is 
essentially the size of the largest cluster of connected ribbons. Measured in units 
of the effective elementary localization length, A = x a , £ diverges as 

I = J_ ~ \AB\-" , (99) 
A xa 

as relative to this length, the problem is that of "ordinary" percolation. However, 
measured relative to a fixed length, such as a, one has that 



- ~ \AB\-^ +1) = \AB\ 



a 



(100) 



and Eq. (95) follows. 

Equation (99) implies that tunneling increases the localization length in com- 
parison to classical percolation, i.e., 

e > ■ (ioi) 

However, Zhao and Feng [46] argues that 

f « e P • (102) 

Their argument goes as follows: We are close to, but not at the classical percola- 
tion threshold, and have no infinite cluster yet. The largest cluster of connected 
paths that would emerge at the percolation threshold has therefore been broken 
up into islands of size £ p . Between each island there is a saddle point where 
tunneling is necessary. The system has a width W, and therefore the number of 
such tunneling saddle points is W/£ p . The probability that an electron tunnels 
through all these saddle points along a path crossing the sample is therefore 



P w/ ip = e -w/(i pX ) = e -w/t ; (103) 



where we have used Eq. (96). The right hand expression defines the localization 
length, and we find 

z = x z p ~\ab\-w> . (io4) 

We have here used Eqs. (92) and (97). From Eq. (97), we see that a tunnel length 
scale may be defined, 



k 



he\AB\ 



(105) 



13 The percolating cluster has the topology of a one-dimensional chain with "blobs" inter- 
spersed among the so-called "cutting bonds." The cutting bonds are those that cuts off perco- 
lation if removed. The blobs are clusters of connected bonds. Coniglio [47] has shown that the 
fractal dimension of the cutting bonds is 1/v at the percolation threshold. 
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and Eq. (104) may be written 

( = (>({)'■ (106) 

Since tunneling has been assumed to be the dominating mechanism for transfer 
across the saddle point, one must have 

k > I ■ (107) 

Combination of this inequality with Eq. (106) leads to the inequality (102). It 
should be noted here that as |A£?| — > 0, this argument implies a crossover to 
classical percolation behavior when l t becomes of the order of /. More detailed 
analytical arguments [48] have been constructed that support the rescaling ar- 
gument in one dimension, but not in two dimensions, which is our concern here. 
On the other hand, no crossover as implied by the argument basic to Eq. (106) 
is found in numerical work. 

In the above arguments, quantum interference has been left out. It would 
lead to an Anderson-type localization [1]. There are no arguments on the level of 
the Mil'nikov-Sokolov that predicts a value for the localization length exponent v 
based on interference. In order to investigate the implications of this mechanism, 
numerical studies are required. We therefore now turn to the numerical modelling 
of the integer quantum Hall problem. 

5. Numerical Models 

In this Section, we turn to the numerical studies of the integer quantum 
Hall effect. Since Mil'nikov-Sokolov-type argumentation is controversial, and no 
arguments taking the full quantum mechanical formulation into account that are 
capable of pinning down the value of u, it seems that the computer will continue 
to be the principal tool for studying this problem. The use of numerical methods 
has, indeed, been quite successful in this field. Not only have results found in 
experiments been well reproduced on the computer, but also computer results 
have had predictive power. For example, the Chalker-Coddington model [42], 
which we will be studying in a moment, predicted in 1988 the valued z> = 2.5±0.5, 
while the experimental result v = 2.3 ±0.1 of Koch et al. [29], was reported three 
years later. 

5.1 The Chalker-Coddington Model 

We will in this review concentrate on network models based on the one in- 
troduced by Chalker and Coddington [42] , as it seems that these models capture 
very well the critical aspects of the integer quantum Hall effect. Other models 

14 Reanalysing the original data of Chalker and Coddington [42] using a different parame- 
trization and including additional data, Huckcstcin [39] found v = 2.35 ± 0.03. 
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Figure 7: Network used in the Chalker-Coddington model. It has length L in the 
horizontal direction and length M in the vertical direction. It is periodic in the 
vertical direction. The amplitudes associated with nodes i and j are shown. 

may be found in Huckestein's review [13]. 

The Chalker-Coddington model replaces the random network of paths in the 
disordered energy landscape by a regular square lattice with a lattice constant 
equal to the average distance between saddle points in the disordered energy 
landscape. The nodes of the lattice represent saddle points, and the bonds the 
paths between the saddle points. Each bond can only be traversed in one direc- 
tion. The possible directions are chosen so that the network forms a grid of loops 
of alternating handedness, see Fig. 7. The wave function itself is represented by 
a complex variable z. At a given saddle point, Z\ and z 2 are ingoing amplitudes, 
and Z3 and z^ outgoing amplitudes as shown in Fig. 8. The scattering matrix is 
conventionally written^] [42] 

(z 1 \_(e i ^ \ /cosh7 sinh7\ /e^ 3 \ f z 4 \ . . 

Us/ V e^V Uinh7 cosl^J V e^ 4 ) \z 2 ) ' ^ ' 

where <pi, <p 2 , 4>3 an d 04, represent the phase changes of the wave function along 
the bonds. These are chosen at random. This randomness represents the varia- 
tion in path length from loop to loop. The control parameter 7 determines the 
scattering properties at the saddle point. When 7 — » 0, the transmission proba- 
bilities are Pi^4 = P 2 ^3 = 1 and Pi^3 = P 2 ^4 = 0. Likewise, when 7 — > 00, we 

15 Traditionally, the scattering matrix is written in the form (108). It could equally well 
have been written in a form that relates the outgoing channels (23, Z4) to the ingoing channels 
(21,22)- 
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a b 

Figure 8: Scattering at a given node in the Chalker-Coddington model: a) shows 
the case when 7 < 7 C while b) shows the case 7 > j c . The dashed curves denote 
the case 7 = j c . 



have that P i _ +3 = P 2 _^ 4 = 1 and P 4 ^ 4 = P 2 _ 3 = 0. When 7 = 7c = ln(l + a/2), 
one finds Pi_> 4 = Pi_> 3 = P 2 ^ 4 = P 2 ^ 3 = 1/2. When 7 < 7 C , Pi^ 3 = P 2 ^ 4 = 
1/(1 + sinh~ 2 7) are due to tunneling, while Pi^ 4 = P 2 _> 3 = 1 — Pl-^ are direct 
channels. When 7 > 7 C , Pi^ 4 = P 2 _> 3 = 1/(1 + sinh +2 7) are due to tunneling 
while the two remaining channels are direct. This is illustrated in Fig. 8. Writing 
P as in Eq. (96), derived in the energy landscape model, we may combine all 
these results in the expression 

^ 2|lnsinh7| ^ 

While Chalker and Coddington [42] kept 7 equal for all nodes, with disorder 
in the phases (pi only, Lee et al. [43] extended the original model by writing 

7 = 7ce"- w , (HO) 

thereby introducing disorder also in the saddle points in addition to the one 
already present. Here \i is a fixed parameter, equal for all nodes in the network, 
while v is chosen randomly on the interval [—w/2, +w/2], where w is a parameter 
determining the width of the distribution. As we will see in Sec. 5.2, the Chalker- 
Coddington model crosses over to classical behavior in the limit w — > 00. In the 
case when \/i — v\ 1, we may relate \/i — v\ oc \B — B c \ of the energy landscape 
model. By using Eq. (98) combined with Eqs. (109) and (110), we find that the 
tunnel length A behaves as 

A = X aK-^—, (111) 
l/i — v\ 
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with a the lattice constant. This expression should be compared to the corre- 
sponding expression, Eq. (97), derived in the smooth-potential model. Thus, the 
renormalization argument presented in Sec. 4.2, predicts the same critical behav- 
ior in the Chalker-Coddinton model as in the energy landscape model. However, 
we reiterate that a localization length exponent v = 2.35 ± 0.03 [39] is found in 
the Chalker-Coddington model without nodal disorder. This has been investi- 
gated further, with the conclusion that the Chalker-Coddington model does not 
support the Mil'nikov and Sokolov argument [49]. 

The way one determines v in the Chalker-Coddington model is based on a 
transfer matrix approach and the assumtion of one-parameter scaling [50]. The 
procedure is as follows. The lattice employed is as shown in Fig. 7, but with 
horizontal length N and vertical length M. We assume N ^> M — typically M is 
of the order of 10 2 and N of the order 10 5 . Calculate the trace of the transmission 
probability from layer one to layer N, tr|T| 2 . Average this quantity geometrically 
over the different samples, i.e., calculate exp[(lntr|T| 2 )]. The reason for doing the 
geometric average is that we want the typical transmission probability. This is 
not accessible with the usual arithmetic average, (tr|T| 2 ), which is dominated by 
the rare events where the sample actually percolates in the classical sense from 
edge to edge (in the M direction), and in which case tr|T| 2 is equal to one. In 
taking the logarithm, these rare events no longer play a role. We may now define 
an effective quantum localization length, £m by 

^- = " i im 4 (lntr|T| 2 ). (112) 

It is effective in the sense that it depends on M. The asymptotic localization 
length, which is the one we want to determine, is then 

e = £m . (113) 

At this point, the one-parameter scaling assumption comes into play. It assumes 
the following functional relationship between £ and £m with M and /x, defined in 
Eq. (110), as parameters: 

T^-F(M] . (114) 
M \ M j y 1 

By plotting for several values of AB and M, data collapse is obtained when 
the correct functional dependence of £ on AB has been found — which is a power 
law with slope —v. 

We shall not discuss the basis for the above procedure here, but instead turn 
to a different way of extracting the localization length exponent from the Chalker- 
Coddington model. One recently employed consists in treating the network as 
an open system [51]. By open, we mean the following: Rather than calculating 
transmission and reflection coefficients from the elementary scattering matrices 
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of Eq. (108), we "hook" the network up to a source and a sink. That is, we 
specify the value of z on all bonds on the edges leading into the network, while 
leaving the z values of the exiting bond unspecified. This corresponds to driving 
a current through the network, and z must now be interpreted as current am- 
plitudes rather than wave functions, as it is no longer normalized. In fact, this 
models the dirty Hall bar with reservoirs shown in Fig. 6, although the boundary 
conditions in the transverse direction are different, v can then be determined as is 
done in connection with the classical random resistor network at the percolation 
threshold.^ 

When the system is considered open, the various link amplitudes are found 
by solving a set of linear equations, 

Az = b , (115) 

where B contains the z values that correspond to specified bonds, and the matrix 
A contains the detailed scattering conditions for each node, based on Eq. (108). 
We fix all currents entering the network shown in Fig. 7 from the left to unity, 
while no currents enter the network from the right. By solving Eq. (115), we 
determine all ^-values in the network, and from these the total current passing 
through the system, namely, 

i= E (kil a -kil 2 )= E M 2 . (lie) 

i £ first column i £ last column 

The first of the two lower indices on the bond variables z refers to the scattering 
channel as shown in Fig. 7, while the second index is the node address. 

The total current entering the network from the left, // is proportional to the 
chemical potential of the reservoir hooked up to the network on that side, while 
the chemical potential on the right is zero, as no current enters from this side. 
Thus, the potential difference across the network is simply V = h/e 2 Ii. We may 
now define a longitudinal conductance G through the expression 

I = VG. (117) 

As is evident from Fig. 7, the nodes belonging to every second row have their 
scattering labels rotated by 90° in the clockwise direction. As a result of this 

16 Suppose the bonds of a standard bond percolation problem are electrical resistors, all 
equal. What are the electrical properties of this network? First of all, all the geometrical 
and topological properties of networks at the percolation threshold are still in place; this is a 
transport problem on top of the original percolation problem. In particular v = 4/3 in two 
dimensions. In addition, one finds that global quantities such as the conductance G scale with 
system size L as 

where t — 1.300 [52]. Its exact value is not known. Furthermore, the current distribution is 
multifractal. We will explain this concept below. See Ref. [53] for a review. 
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rotation, the connection between bond labels and forward/backward scattering 
depends on which row the nodes belong to. For the odd-numbered rows, scat- 
tering z\ —>■ Z3 represents forward scattering, whereas for the even rows, it is 
represented by z-i — > z$. The scattering amplitude for channel z\ — > z% increases 
with increasing 7, whereas the scattering amplitude for channel Z2 — > z% increases 
with decreasing 7. 

Rather than assigning random values for 7 to the nodes in the network, we 
have used two bimodal distributions, one for the odd rows, 

Vo{l) = (1/2 + n)8in - 713) + (1/2 - n)8{>) - 723) , (118) 

and one for the even rows, 

p e ( 7 ) = (1/2 - ^(7 - 7l3 ) + (1/2 + /i)5( 7 - 723 ) . (119) 

These two distributions contain three parameters: 713, 723 and /1. We relate 713 
and 723 through the relation 

tanh 7l3 = — I = t , (120) 

cosh 723 

where to is a transmission amplitude. With the relation (120), a node belonging 
to an even row with a 7 = 723 will transmit in exactly the same manner as a node 
belonging to an odd row with 7 = 713. When |to| — * 1, tunneling becomes less and 
less likely and the model becomes classical, as will be discussed in Sec. 5.2. On 
the other hand, in the |io| — ► 1/v^ limit we find that j 13 = 723 = 7 C = ln(l + \/2). 

With a large negative value of control parameter n, the system will predomi- 
nantly consist of elementary loops with current running clockwise, whereas with 
a large positive value of fi, the current will be running counterclockwise in loops 
shifted by one lattice constant in the (11) direction. When /i = 0, there is no 
bias and the system is critical. 

In our simulations, we set to — 0.9, and let fi vary between -0.5 and 0.5. In 
Fig. 9 we show the current I as a function of \x for square system sizes L = 12 to 
92. The data were geometrically averaged over 14650 samples for L = 12 to 300 
samples for L = 92. 

From Fig. 9 we see that the larger the system size, the sharper the peak. The 
system becomes critical when L/£(A/x) ~ 1, hence we find that the half width, 
A/i, scales as 

Afi oc L- l/i> . (121) 

By plotting the half width A/i versus system size in a log-log plot, Fig. 10, we 
obtain the localization length exponent v = 2.35 ± 0.07.^] 

17 We have also performed similar simulations in the limit when to = 1. As these simulations 
are much less computer intensive, we have averaged up to 100 000 samples of systems ranging 
in size between L = 10 and L = 960, obtaining v — 1.333, in accordance with the expected 
value 4/3 — see Sec. 5.2. 
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Figure 9: The current passing through the system as a function of chemical 
potential, for various system sizes. 

So far, we have only discussed one single critical exponent associated with the 
localization-delocalization exponent in the quantum Hall effect, namely u, the 
localization length exponent. Recent analysis of the Chalker-Coddington model 
by Klesse and Metzler [54], which is based on earlier work of Pook and JanBen 
[55], shows that the wave function is multifractal in this model. In terms of an 
open system, where z is to be interpreted as a current amplitude, the problem 
is operationally analogous to the random resistor network at the percolation 
threshold, where it also has been established that the current distribution is 
multifractal [56,57]. 

Multifractality means the following: At the critical point, the logarithmically 
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binned histogram, H, of currents z show the following scaling form: 

H(z,L)~L f ^, (122) 

where 

z ~ L~ a , (123) 

and L is the lattice size (assuming a square lattice). The two functions / and a 
are independent of L. If f(a) is a more complicated function of a than of linear 
form a + 6a, then the current distribution is multifractal. This signifies that 
there is a continuum of singularities governing the distribution — singularities in 
the sense of non-trivial power laws, and not just a small number of exponents, 
which is usually the case in connection with critical phenomena. The geometrical 
scaling properties of percolation, for example, are completely described by two 
exponents, while the current distribution is multifractal.0 

One can study a distribution by building a histrogram. Equivalently, one can 
explore it by measuring its moments. In our case, we define the nth current 
moment as 

M n = (^|z| 2n ) ~ L yM . (124) 

i 

They all turn out to scale as power laws. The tell-tale sign of multifractality is 
that 

y{n)^y(0) + cn, (125) 

where c is a constant. The fundamental relation between /, a, y and n is of the 
Legendre type, 

y(n) = f — na . (126) 

This interpretation and equation form the heart of the celebrated f — a formalism 
[58]. 

We will not pursue this line any further here, except for pointing out that this 
is today an active field of research. 

5.2 Classical Limit of the Chalker-Coddington Model 

The Chalker-Coddington model has an interesting classical limit. For finite /x 
and v — > oo, defined in Eq. (110), the scattering channels are z\ — > z± and Z2 — > Z3, 
and when v — > — 00 they are z\ — > Z3 and Z2 — > 24. In this extreme case, there is 
no contact across the saddle points. However, derealization is still possible. The 
localization length exponent in the extreme case when w — > 00 has been measured 
by Lee et al. [43] to be 1.29 ±0.05 in agreement with classical percolation theory 
(predicting 4/3). However, Bratberg et al. [51] have noted that the system in this 

18 To clarify Eq. (122) and (123) somewhat, think of a (i-dimensional resistor network without 
percolation disorder: All bonds are present. Then, a — 2 — d and f = d. The / vs. a curve is 
in this case just a point. 
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Figure 10: Log-log plot of the half width of the curves shown in Fig. 9. The 
straight line is the result of a least squares fit giving v = 2.35 ± 0.07. 

limit is not equivalent to the percolation picture presented in Sec. 4, but to the 
system shown in Fig. 11. This is a random tiling of the two elements shown at the 
bottom of the figure. The curves that the markings on the tilings — known as 
Truchet tiles, see e.g. Gale [59] — produce are known as Dragon curves, see Wells 
[60]. The critical properties of this tiling was studied some time ago by Roux et 
al. [61]. Their conclusion was that it is in the universality class of smart kinetic 
walks (SKW). These may be mapped on the external perimeters of percolation 
clusters in two dimensions. As a consequence, their correlation length exponent 
is the same as in percolation: vskw — v — 4/3. Thus, the Chalker-Coddington 
model with its good exponent v = 2.35 ± 0.03, does not simulate the percolation 
picture of the localization-delocalization transition in the integer quantum Hall 
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Figure 11: The Chalker-Coddingon model in the limit w — * oo generates "dragon 
curves" by randomly tiling the plane with the two Truchet tiles shown. 

system in the "classical" limit, w — > oo, where w is the width of the distribution 
from which v is chosen. The criticality of the model in this limit comes from the 
loop size distribution, and not from joining of overlapping loops at saddle points, 
as there are no saddle point crossings. 

6. Conclusion 

We summarize in a few sentences the picture of the integer quantum Hall effect 
that we have presented here. Non-interacting electrons move two-dimensionally 
in a strong perpendicular magnetic field and in a disordered potential. Clas- 
sically, the perpendicular magnetic field causes the electrons to move in circles 
about the so-called "guiding centers." Quantum mechanically, this circular mo- 
tion is quantized into Landau levels. We assume the disorder potential to be 
slowly varying on the scale of the radii of the Landau levels, so that there are no 
transitions between these caused by the potential. The Schrodinger Eq. (45) then 
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simplifies to Eq. (49), which shows that the energetically allowed regions where 
the electron may move are the equipotential curves of the disordered potential. 
The wave function is localized to bands of width of the order of the radius of the 
Landau level, and follows the equipotential curves of the potential. By adjusting 
the perpendicular magnetic field, these curves move up or down in the disordered 
potential landscape. At a particular critical value of the magnetic field, there 
will be a continuous path of equipotential curves that span the system. The 
electrons in this Landau level are then delocalized, and there is longitudinal con- 
ductance. For other values of the magnetic field, they are not. This mechanism 
is repeated for every filled Landau level up to the Fermi level. The localization- 
delocalization transition is driven by coalescence of disjoint equipotential curves, 
and strongly resembles classical percolation. However, tunneling at the saddle 
points and interference effects change the universality class of the transition from 
that of percolation. There are semi-classical arguments that reproduce the der- 
ealization length exponent seen numerically and experimentally, which essentially 
describe the localization-delocalization transition as a "dressed" percolation pro- 
cess. However, these arguments have serious weaknesses, and it seems likely that 
the transition is a pure quantum one. 

The number of experiments that directly test the scaling properties of the 
transition are few and to a certain degree not consistent with each other. Why 
do for instance Dolgopolov et al. [35] find a localization exponent v ps 1 which 
is roughly consistent with the classical value v = 4/3, rather than the quantum 
value w 2.3? The argument used to derive the exponent v ~ 1 in the work 
of Dolgopolov et al. [34] is equivalent to the gradient percolation argument of 
Sapoval et al. [62], which is in fact incompatible with the slowly- varying landscape 
model presented in these lectures. 

This is clearly not a closed subject. 
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